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Over the last several years the interest in the class of 
magnetic Bp stars has spurred the accumulation of datasets 
at many wavelengths over the stars' rotational/magnetic 
cycles. These datasets support the view that the circum- 
stellar environment is the site of a balance of physical pro- 
cesses that channels the wind toward the magnetic plane. If 
the magnetic field is strong, that is if the circumstellar mag- 
netic energy density exceeds the wind energy density, the 
star's wind is channeled into the magnetic plane where a 
shock is formed and the cooled shock settles to form a stable 
disk. The magnetic axis is typically inclined with respect to 
the rotational axis. Then, as different portions of the disk 
wobble in front of the star during the rotational cycle, they 
produce alternate absorption and emission components in 
the Ha, UV resonance lines, and low-excitation metallic 
lines (e.g., Bohlender et al. 1987, Shore 1987, Bolton 1994). 
To complicate the picture, it appears that slow infall of disk 
matter is responsible for producing inhomogeneous metal- 
poor patches or a belt along the magnetic equator (e.g., 
Khokhlova et al. 2000, Groote 2003). Also, although there 
remains a theoretical difficulty with the fractionation hy- 
pothesis for helium (Krticka et al. 2006), enough neutral 
helium atoms in fact seem likely to separate from the polar 
wind, return to the surface, and accumulate there as helium 
caps. 

Recent theoretical refinements in the picture (e.g., 
Preuss et al. 2004, Townsend & Owocki 2005; hereafter 
TO05) suggest that the wind will settle onto low equipoten- 
tial surfaces determined by the balance of radiative, grav- 
itational, magnetic, and centrifugal forces. These surfaces 
reside primarily, though not exclusively, in the disk plane. 
Because the disk is likely to have a nonaxisymmetric den- 
sity distribution, we will refer to the accumulations near the 
plane responsible for time- variable spectral features as disk 
segments. TO05 have presented a "rigidly rotating magne- 
tospheric" ah initio model to explain physical properties of 
these segments for the case of ct Ori E and it can be applied 
to the cases of other Bp stars with arbitrary magnetic incli- 
nations and viewed from any rotational inclinations. TO05 
find that the disk plane is not precisely perpendicular to the 
magnetic polar axis for intermediate magnetic obliquities, 
including those they obtained for crOriE. 

The spectroscopic anomalies of cr Ori E (HD 37479, 
B2Vp; Lesh 1968) were first noticed by Berger (1956) and 
discussed further by Osmer and Peterson (1972). This star 
is the prototype and most studied of the helium-strong sub- 
class of Bp stars. Walborn (1974) discovered Ha emissions 
in this star's spectrum that varied on a timescale of hours. 
This variability was subsequently established to be asso- 
ciated with the rotation of a frozen-in disk, as described 
above. Hesser et al. (1977) determined a period of 1.19081 
days based on several nights of photometric observations 
distributed over a few years. Subsequently, various studies 
(Landstreet & Borra 1978, Bohlender et al. 1987) deter- 
mined that the magnetic variations (suggestive of an ap- 
proximately dipolar field with polar field strength Bp « 
10 kG) are consistent with the photometric period. Various 
periods over the range 1.19081 — 1.19084 days have since 
been determined by several authors from analyses of vari- 
ations of photospheric and circumstellar Ha line emissions 
(Hesser, Moreno, & Ugarte 1977, Reiners et al. 2000, here- 
after ROO; Townsend, Owocki, & Groote 2005; hereafter 
TOG05). 



Groote & Hunger (1976, 1977 hereafter GH76 and 
GH77) discussed the variations of high-level Balmer lines 
based on only photographic spectra. The GH76 data in- 
dicated that two disk occultations of the star, about 0.4 
cycles apart, cause high-level Balmer line absorptions visi- 
ble at these times. At the same times, the intervening disk 
material absorbs flux just shortward of the Balmer edge. 
Although the analysis of line absorptions do not generally 
allow us to infer the extent and total volume of such a 
corotating structure, these absorptions do yield consider- 
able information about its three dimensional geometrical 
character. Smith & Groote (2001, hereafter SGOl) com- 
bined strength and shape information of many ultravio- 
let metallic lines obtained by the International Ultraviolet 
Explorer around the rotation cycle to exploit this poten- 
tial. Their study provides some of the first estimates of the 
physical parameters of the disk gas, including disk temper- 
ature, column density, and areal coverage of the star. The 
limitations of this study were the low signal-to-noise ratio 
and the paucity of phase sampling of the archival database. 
What has been needed is a study that includes complete 
rotational phase coverage of absorption lines sensitive to 
electron density. To respond to this need, we have obtained 
complete phase coverage of the hydrogen lines from H9 
(HC) to the Balmer edge over the entire rotation period 
of CT OriE. Such datasets have the potential to better eluci- 
date areal coverage and density structure of the absorbing 
disk segments than has previously been possible. 



1. Observations and modeling methodology 

1.1. Observations 

In an attempt to search for variations in high-level 
Balmer lines from the variable occultation of a Ori E by 
its circumstellar disk, we obtained 114 spectra with the 
1.85 m Plaskett telescope of the Dominion Astrophysical 
Observatory (DAO) on seven nights between 2007 January 
25 and February 28. We used the f/5 Cassegrain spectro- 
graph with an 1800 1/mm grating blazed at 5000A and 
the SITc-2 CCD. Despite the sometimes mediocre seeing 
at the DAO, we used a narrow 1 arcsec slit rather than an 
available image sheer in order to maximize the resolution 
and to permit background sky subtraction. We note, how- 
ever, that observing conditions and seeing were very good 
on each of our observing nights. Exposure times were 15 
minutes and to minimize the effects of any possible flex- 
ure in the spectrograph Fe/Ar comparison arc reference 
spectra were obtained approximately every hour. The data 
were processed in a standard fashion with IRAFQ The re- 
sulting spectra have a dispersion of 10 Amm~^, a reso- 
lution of about 16400, and cover the wavelength region 
AA3645— 3863. Typical signal-to-noise ratios, as measured 
in a flat region of the spectrum near A3850, are about 150 
per pixel. Table 1 provides a journal of observation numbers 
(the DAO 1.8-m odometer numbers for 2007 spectra), ob- 
serving times expressed as Heliocentric Julian Day (HJD) 
and rotational phases. 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy (AURA), Inc., under 
cooperative agreement with the National Science Foundation. 



2 



Table 1. crOriE observing log including DAO 1.8-m spectrograph 2007 odometer numbers. Heliocentric Julian Dates - 
2450000.0 and rotation phases are given in columns 2 and 3. The phases are derived using a formal period of 1.190841 
days and the 1976 zeropoint of HJD 2442778.819 used by Hesser et al. (1977) - see text. 
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1.2. Spectral synthesis computations 

To perform the analysis of the high-level Balmer line and 
continuum spectra, we utilized a pair of programs written 
by I. Hubeny and T. Lanz with both LTE (Kurucz 1993) 
and non-LTE models in the BSTAR2006 scries (Lanz & 
Hubeny 2007). The first of these, synspec, is a photospheric 
line synthesis program that computes the line and contin- 
uum photospheric spectrum of a star (Hubeny, Lanz, & 
Jeffery 1994). We also made use of synspec^s capability of 
convolving the photospheric lines with functions approx- 
imating the instrumental and rotational broadening. We 
computed the synthetic spectra in steps of O.OlA. 

The goal of our analysis was to compute effects of the 
high-level Balmer lines and Balmer continuum due to the 
intervening disk during each of the two occultation events 
that occur during (jOriE's rotation cycle. Since it is al- 
ready known from the SGOl study that the disk is cooler 
than the photosphere, these effects will only include absorp- 
tions. To simulate these absorptions we used the radiative 
transfer program circus (Hubeny & Heap 1996). This pro- 
gram computes strengths of absorption or emission com- 
ponents of lines in a circumstellar medium from user-input 
quantities such as disk temperature, column density, arcal 
coverage factor (the portion of the disk that passes in front 
of the star), and inicroturbulent velocity. Circus combines 
this contribution with the st/nspec-computed spectrum of 
the photosphere. In the solution of the radiative transfer 
in its "LTE mode," circus calculates the line absorption 
reemission along a line of sight according to an input disk 



temperature, Tdisk- The physical rationale for choosing this 
convenient mode is that the formation of the Balmer line 
flux in the disk is believed to be due to recombination, and 
the atomic levels of the high-level lines and continuum may 
be plausibly assumed to be close to equilibrium with the 
gas kinetic temperature. As a reference point, we first con- 
sidered temperatures near 12 750 K, the value SGOl found 
for Tdisk from an analysis of low-excitation UV metallic 
lines. As detailed below, at one point in our analysis we 
also used the circus "scattering" (no reemission) approxi- 
mation to gauge the effects of non-LTE in the disk forma- 
tion of Balmer line cores. In addition, we point out that 
absorptions from the corotating disk are formed at these 
same Doppler velocity as for the background star, i.e., no 
differential velocities were imposed along the line of sight. 
Although circus is capable of computing the spectrum for as 
many as three circumstellar cloud components, our initial 
assessment of the equivalent width (EW) variations sug- 
gested that the absorbing disk segments have a complicated 
density stratification perpendicular to the disk plane. Since 
these details can not yet be explained by ab initio models, 
it was sufficient to describe the absorptions in terms of a 
single-component cloud model in our circus analysis. 

2. Description of disk-induced absorptions 

2.1. Equivalent width variations over the cycle 

To survey the absorptions of the high-level Balmer scries 
in our spectra, we chose the members H9, H12, and H14 



3 



for EW measurement. To these members we added the 
HeA3819 hne, as it showed important differences with re- 
spect to the Balmer lines. We selected the Balmer H9 line 
for study because of its proximity to the He I line and there- 
fore the convenience of exploiting a nearly common local 
continuum level that would allow reliable measurements 
of line core and wing variations. The HI 2 and H14 lines 
were chosen arbitrarily because they measure effects of yet 
higher atomic levels (such as diminished lino opacity), and 
they are largely redundant to one another. This advantage 
allows one to check the consistency of the EW variations. 

For our line analysis we first rectified each spectrum by 
identifying maximum fluxes along the spectrum and fltting 
consistent orthogonal polynomials through these points us- 
ing the IRAF task continuum. "Raw" EWs were measured 
from fluxes within ±2.5A of the line centers; for H14 this 
window was relaxed to ±2.4A. These windows were chosen 
because they delimit the absorption from the disk occulta- 
tion. For the initial measures these core depths were refer- 
enced to a pair of quasi-continuum points intermediate be- 
tween the line to be measured and the nearest hydrogen and 
helium lines. One of these points falls between the closely 
occurring Hcl A3819 and H9 lines. This circumstance al- 
lowed us to check for residual errors in local continuum 
placement. Such errors would produce similarly low or high 
EW measures in both these lines at the same phase, rather 
than the opposing variations in these two linos caused by 
deviations from the mean of the surface He/H value as the 
star rotates. Because our analysis is based on absorption 
excesses of the line cores, our measurements can be tied 
directly to the computed absorptions in our disk models. 

We noticed that our raw EWs were significantly af- 
fected by an approximately single-waved variation in the 
line wings that is well correlated with the variation of 
the He I line measures and hence instantaneous integrated 
He/H abundances measured by ROO. We corrected for this 
effect by first measuring the variable fluxes in their wings at 
windows between ±2.5A and ±5A from line center. These 
fluxes are not affected by the disk. The wing EW compo- 
nents of the H and He I lines vary in perfect anti-correlation. 
We chose an arbitrary reference phase of </) = 0.35, which 
corresponds to the middle of an extended low plateau of 
the strength of helium lines (ROO), signifying a near- normal 
He/H abundance for the region of the star visible during 
this phase. We then corrected the core EWs using an empir- 
ical relation derived from the slowly varying wing and core 
variations for the H9 and He I lines throughout the cycle 
(excepting the brief occultation maxima discussed below). 
Noticing that this correction resulted in an equality of the 
two H9 lino EW minima, we adopted similar relations for 
the H12 and H14 line EWs that forced equal EWs in their 
minima as well. Figure 1 exhibits the corrected variations 
through the rotation cycle for these four lines. Different 
symbols are used for different nights' observations to best 
exhibit the (generally negligible) night to night differences 
in our measurements at thoses phases observed on multiple 
nights. 

The panels exhibit several interesting details in the 
Balmer and He I line absorption curves of crOriE. First, 
the EW curves exhibit two maxima at (j) ^ 0.0 and 0.4. 
Distinct minima are also apparent at </> = 0.85-0.90 and 
0.55. We will refer to these as the primary and secondary 
absorption maxima and minima, respectively. The excep- 
tion to this statement is that the primary maximum is ab- 



sent in the He I line curve and the secondary maximum is 
only weakly present. The existence of the brief minima was 
surprising because they have not been noticed in previous 

EW studies in the UV and optical (e.g., ROO, SGOl). As 
our circus analysis will also show, all of the lines are opti- 
cally thick, but the H14 line is significantly less thick than 
H9. Additionally, one notices differences in the approach 
of these curves to their maxima and minima. The edges of 
some of these plateaus arc sharp, and this is likely to be a 
manifestation of the line opacities changing from optically 
thin to thick regimes. 

The primary maximum differs in morphology from the 
secondary feature, being more short-lived and triangular. 
Observation ^800 coincides with the maximum EW for all 
the H lines. Following, Hesser et al. (1977), we assign </> 
= 0.0, defining the line flux minimum with the midpoint 
of this observation #800, HJD 2454159.687, as the epoch- 
defining flux minimum. For completeness, we note that if 
one adopts Hesser et al.'s 1976 epoch of HJD 2442778.819, 
uses the cycle counts of ROO, and assumes that the period 
has been constant over this 10,000 cycle interval, a period of 
1.190841 days would result. However, we hesitate to claim 
this formally as the star's true period (or to assign error 
bars) because the apparent phase difference in photometric 
minimum and line strength maximum of 0.02 cycles sug- 
gests that caution should be exercised in combining heter- 
gencous data types - that is our (p = corresponds to 0.02 
for the Stromgren filter u minimum phase of Hesser et al., 
according to the GH77 phases. If correct, this would mod- 
ify the value of the period to 1.190839 days. We also point 
out that the most current value of a wholly photometrically 
derived period, 1.190832 days (Townsend 2007) is only 3(t 
from the above value, according to the errors we would per- 
haps simplistically assign to our own measurements. 

In addition to the maxima and minima, minor brief 
changes in the absorption take place, such as that occurring 
at (f) = 0.80 in the strong lines H9 and He I A3819. More 
conspicuously, the flat plateaus at = 0.7—0.8 in the EW 
curves of these lines are not shared in the weaker HI 2 and 
H14 curves, suggesting a thinning of absorbing circumstel- 
lar material that is first betrayed in the least optically thick 
lines. Brief minor absorption bumps can be seen in these 
curves near <p = 0.20 and 0.93. Altogether, there is sub- 
stantial evidence for the existence of static, high-latitude 
circumstellar matter. The fact that these absorptions are 
present to some extent in He I argues as well that this mat- 
ter is heated more than, for example, the disk segment re- 
sponsible for the primary absorption. The excitation of this 
matter is still influenced by a largely undiluted UV radia- 
tion field from the nearby star. 

We end this section by discussing whether the EW vari- 
ations in Fig. 1 could be due to the passage across the disk 
of surface inhomogeneities, which could influence the line 
strengths either through changing the number of absorbers 
in the line or indirectly through their influence on the struc- 
ture of the photosphere. To address this hypothetical prob- 
lem, we stress that the wings of these lines do respond to 
such changes consistently with EW variations of the He I 
lines including the A4713 line depicted by ROO (e.g. dur- 
ing the phase interval (j) = 0.65—0.80) and do not deviate 
from the anti-correlation established from this interval at 
intervals when the HI and Hel line cores undergo sudden 
changes. Second, the occultation maxima found in Fig. 1 
develop in very short intervals, contrary to the slower vari- 
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Fig. 1. Equivalent width curves plotted as a function of rotational phase using the ephemeris in Table 1. The EWs are 

defined from local "continuum" points located between the hydrogen lins. Represented are panel a: H9, panel b: H12, 
panel c: H14, and panel d: Hel A3819 lines observed during 2007. Symbols are: pluses (25 Jan), stars (29 Jan), diamonds 
and dots repeated a cycle earlier (30 Jan), triangles (31 Jan), squares (1 Feb), crosses (2 Feb), and diamonds (28 Feb; at 
^ = 0, corresponding to observation #800 for all lines). 



ations of at least 0.15 cycles for which the Hel and metal 
line strengths evolve (ROO, Groote 2007). 

A third reason for doubting that inhomogeneities influ- 
ence our line wing measurement derives from results using 
specialized model atmospheres results appropriate to Bp 
stars with He and metallic patches. In particular, Krticka 
et al. (2007) have studied the effects of chemical inhomo- 
geneities on the surface of the B2p star HD 37776, another 
B2p star with a very strong magnetic field. Both the abun- 
dance variations across the surface and the magnetic field 
strength arc at least 10 times greater than those in a OriE. 
These authors find that photometric variations of this star 
are likely due to the effects of bound-free absorption edges 
produced by abundance patches of elements like silicon (the 
iron abundance is depleted). They show that the temper- 
atures of superficial regions of the photosphere should be 
elevated by this added opacity. Because this is the gen- 
eral region where most line cores and bound-free edges are 
formed, the effect is to cause a weakening of the line cores 
arising from "passive" species, in this case helium and hy- 
drogen atoms. The metals and helium variations are anti- 
correlated because for both a Ori E and HD 37776 the met- 
als are less depleted in belts distributed on the magnetic 
equator. Thus, the passages of patchy regions of the star's 
surface, irrespective of an intervening disk, produce contin- 



uum light variations. Returning again to our spectroscopic 
data for a Ori E, and utilizing the Krticka et al. simulations, 
the cores of hydrogen and helium lines should weaken dur- 
ing these phases - yet they do not. We note in particular 
that this fact shows that the slight enhancement of the He I 
A3819 line in Fig. 1 cannot be attributed to surface inho- 
mogeneities on the star. 

2.2. Line profile variations 

To represent the line profile variations of the two occul- 
tations, we compared individual or pairs of line profiles 
differenced against average profiles of the two minima. 
Specifically, for the first stretch of absorption minimum cen- 
tered at (j) — 0.87 the mean was computed from observations 
#305—311, and for the second minimum (cj) = 0.55) from 
observations #448—457. Figures 2 and 3 exhibit, at the bot- 
tom of each panel, the migration of absorption lobes across 
the profile for Hel and H9, HI 2 and H13 line, and HI 5 
and higher for both occultation events. At the top of each 
panel we depict the rectified unabsorbed mean spectrum 
and an example of one of the observations during the oc- 
cultation. Close inspection of the difference spectra during 
the second occultation discloses the passage of two lobes 
across the profile. In our depiction the two lobes can be 
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Fig. 2. Montage of mean and differences profiles through the primary occultation event. Bottom spectra depict difference 

profiles with time (downward); top spectra are the mean spectrum of the non-absorption times at (/) — 0.85—0.90 and 
the observation ^802 (cj) = 0.18) for which the central absorption lobe crosses the mean profile. Panels a, b, and c 
represent the spectral region around H9 and A3819, H13 and H14, and H15— H24. Sloped dotted lines trace the migration 
of subfeatures in three lines with time. The difference spectra refer to observations #798 & 800, 801-802, 803, 805-806, 
and 807-808. These correspond to approximate phases of 0.99, 0.01, 0.03, 0.04, and 0.06, respectively. 



discerned in Fig. 3 by the flat-bottomed difference profiles 
of observations #513-515 (0 = 0.31-0.33), followed by a sec- 
ond flat-bottom in the #523-524 {(f) = 0.38-0.39) profiles. 
The second lobe in the second occultation can be param- 
eterized by a ~50 kms^^ gaussian and is deeper than the 
first one. The first lobe of this occultation, as well as the 
lobe of the first occultation, are slightly narrower and can 
be fit with a ~40 kms~^ gaussian. 

From Fig. 2 we estimate that the midpoint of the first 
occultation, as judged from the transit across the profile of 
the centroid of the excess absorption lobe, occurs for ob- 
servation #802 {(f) = 0.018), or fti0.02 cycles later than the 
observation responsible for the maximum hydrogen EWs in 
Fig. 1. Our error estimate given above was obtained from 
the 0.02 cycle difference between the light minima and line 
profile passage events found by GH76. 

The difference between times of passage of the absorp- 
tion line in the profile and EW maximum for the whole 
line carries over to the second occultation. From Fig. 3 be- 
cause of the finite signal to noise, the breadth and unequal 
depths of the features we cannot establish a reliable phase 
for the mid-point of the occultation. Our best estimate is 
at the beginning of observation #524 = 0.390), which 
would be some 0.37 cycles after the central lobe passage in 



the primary event and identified with observation #802 (0 
= 0.018). Thus, whereas there is a separation of 0.40 cy- 
cles between the EW maxima, the separation between the 
lobe passage events is somewhat smaller. Similar structure 
can be discerned in the semi-resolved double lobes, centered 
at (/) « 0.35 and 0.42 of the secondary occultation in the 
light curves of Hesser et al. (1977) and Oksala & Townsend 
(2007). Perhaps related to the complicated shape of the sec- 
ondary EW maximum is the fact that our best estimate of 
the separation between this and the primary occultation, 
0.40 cycles, is about 0.02 cycles smaller than the value de- 
termined from the Hesser et al. (1977) and the nearly con- 
temporaneous Oksala & Townsend (2007) light curves. It is 
not clear to us whether this minor discrepancy has any con- 
sequences on the details of the disk model. The most con- 
servative assumption is that the opacity (column density) 
of the disk is unequal at opposing positions on either side 
of the central plane. As already suggested from resonance 
line data (GH97, SGOl), a disk "warping" could explain in 
principle both the appearance of a minor lobe within the 
secondary occultation event as well as an apparent unequal 
fall off in column density at small and opposite phases from 
the central maximum absorption. This circumstance would 



6 



3820 3830 3840 3850 




3670 3680 3690 3700 3710 

Wavelength (A) 



Fig. 3. Montage of mean and differences profiles through the secondary occultation event. Pairs of difference spectra with 
time (downward) are depicted at bottom. Spectra at the top are the mean spectrum of the non-absorption times at « 
0.55 and the observation #524 for which the central absorption lobe crosses the mean profile. Difference spectra refer to 
observations # 509-510, 513-515, 516-518, 519-520, 521, 523-524, and 525-526. These correspond to approximate 
phases of .29, .32, .34, .36, .37, and .38, and .40, respectively. 



likely affect optically thin (photometric) and thick (Balmer 
line) measurements differently. 

In our interpretation, this implies asymmetric density 
distributions above and below the magnetic planes; that is, 
this disk segment has partially separated layers. 

In the curves of the primary and secondary occulta- 
tions (Figs. 2c and 3c) the Balmer line series can be dis- 
cerned out to H24 and II23, respectively. (We note that 
differences formed against the mean of all spectra except 
those in the two occultations have a higher signal quality 
and bring out these limiting lines more clearly.) We will 
return to these limiting line members in our later discus- 
sion of electron densities in 331 In the lower panels of Fig. 2 
and 3 we have also indicated the acceleration of the absorp- 
tion lobes through the profile. A cross-correlation analysis 
of the shifts of the line with respect to a fiducial spectrum 
discloses that the average acceleration of these features is 
72 kms~^ hr~^ and 50 kms~^ hr^^ for the primary and sec- 
ondary occultations, respectively. The larger accelerations 
of the features occurring during the primary occultation 
(assuming the disk is rigidly rotating) indicate that it lies 
further from the star than the disk segment producing the 
secondary occultation. As pointed out earlier, the absence 
of the He I line absorption during the primary occultation 



is likewise consistent with this segment lying further from 
the star. 



2.2.1. Balmer line and continuum absorptions 

In preparation for a circus analysis of the physical parame- 
ters of the absorbing disk segments, we have computed ra- 
tioed spectra of the means of the raw (unrectified) spectra. 
These exhibit the largest absorptions during the two occul- 
tations (observation #800 for the primary, and observations 
#516-524 {(j> w 0.36) for the secondary) with respect to the 
absorption minima spectra occurring at nearby phases (ob- 
servations #305—311 and observations #448— 457, or (j) « 
0.86 and 0.56, respectively). These spectra and the ratios 
are depicted in the two panels of Figure 4. These panels 
show the clear absorptions in both the lines and the effec- 
tive Balmer continuum, which starts at A3690— 3700. The 
ratios F(A3650)/F(> A3700) for the occultations come to 
about 0.845 and 0.925. We will use circus to reproduce these 
ratios in the next section. 

It is worth pointing out that the Balmer continuum ab- 
sorption in these particular pairs of specra are unlikely to be 
contaminated significantly by the effect of telluric UV ab- 
sorption through different air masses. For example, in the 
case of the primary occultation, the observation was ob- 
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Fig. 4. The spectrum exhibiting the maximum absorption 
during central occupation and the mean spectrum repre- 
senting no disk absorption referenced in previous figures. 
The ratio of these spectra are depicted at the top. Panels 
a and b represent the primary and secondary occultations, 
respectively. In a, by "M" we denote metallic line features 
used in the text to set a rough disk temperature and also 
the weakened helium line (shown as apparent "emission") 
due to chemical inhomogeneities. The Balmer jump is rep- 
resented by the lower flux shortward of about A3690. 



tained at an hour angle (HA) of 0:49, as compared to the 

nearly equivalent value of a mean HA of 0:47 for the ob- 
servations comprising the iniiiimum spectrum. The respec- 
tive HAs for the secondary event are likewise practically 
identical: 0:33 and 0:36. These differences are less than the 
exposure times of the spectra (15 minutes). 

The hydrogen line absorptions show no noticeable ra- 
dial velocity shift (< ±5 kms~^), and they are confined 
entirely within ±2. 5 A of the line cores. One anomaly in 
the primary occultation spectrum (Fig. 4a) is that the He I 
A3819 actually appears ?/;eaA;er during the occultation. This 
effect reflects the passage of a He-rich spot at phases (j) = 
0.8-0.9 (ROO), from which spectra used to compute the hy- 
drogen line minimum were taken. 



3. Determination of disk parameters 

The primary estimates we are able to make for the disk 
segments causing the hydrogenic absorptions come from 
analysis with the synspec and circus models. We adopted a 
solar H/He abundance in our disk models. However, since 
these are unlikely to be strictly correct (assuming the He/H 
abundance in the disk sustantially reflects the values for the 
material at both magnetic poles), the column densities we 
derive should be understood to refer to hydrogen atoms. 

3.1. Rotational and magnetic axes inclinations 

We have discussed the likelihood already from the discus- 
sion of EW and line profiles that a pair of disk segments 
residing near the intersections of the rotational and mag- 
netic equators occult different regions of the star at times 
separated by approximately one-half of the rotation period. 



Our disk analysis begins with an estimate of the stellar 
properties. We used Kurucz (1993) models with parame- 
ters TefF = 23 000 K and logg = 4 (e.g.. Smith & Groote 
2001). We computed the photospheric line profiles with a 
rotational velocity of 150 kms^^. The foregoing geometric 
properties of the disk are not sensitive to these estimates. 

The geometric constraints begin with a determination 
of the inclination angle i of cr OriE's rotational axis as 
well as the angle between the rotational and magnetic 
axes - the so-called magnetic obliquity, /3. Traditionally 
the former comes from comparing the star's projected ro- 
tational velocity with its radius as inferred from its posi- 
tion on a stellar evolutionary track. Estimates in the liter- 
ature range from i « 90° (Bohlcnder et al. 1987, Short 
& Bolton 1994, TOG05), to i « 39° (ROO). The latter 
value is based on an uncertain distance estimate. Prom de- 
tailed fitting of the light curve of the star to their Rigidly 
Rotating Magnetosphere (RRM) model, TOG05 found that 
/3 + i sd 130°. They adopted i w 75° and /3 « 55°. The in- 
termediate value of /3 is based largely on the unequal spac- 
ings (0.4, 0.6 cycles) of the light minima of the light curve, 
which implies that the center of the stellar magnetic field 
is displaced 0.3R, perpendicular to both the rotational and 
magnetic axis and is also displaced 0.05R, along the mag- 
netic axis. We note also that Groote (2003) gives /? = 67°. 
However, this value is based on an analysis of the distribu- 
tion of He/H on the surface. In particular, it is likely that 
the centroid of the He-rich patch is somewhat displaced 
from one of the two magnetic poles (e.g., Groote & Hunger 
1997). An additional set of constraints can be placed on 
these angles from the shapes of the equivalent profiles of 
the two maxima. As already mentioned, the primary max- 
imum in Fig. 1 has a triangular shape, and the two legs 
of this triangle have unequal durations. This occultation 
appears to be the one for which the segment passes most 
closely to the projccitcid center of the star in the observer's 
line of sight. At (p = 0.4 the opposing segment occults only 
the lowermost limb of the star during mid-passage. We will 
now investigate the shapes and intensities of the occultation 
profiles. 

We have designed and discussed a disk occultation tool 
to quantify the absorption of light and equivalent width 
curves from an optically thick disk or, alternatively, from 
a thick disk surrounded by a translucent periphery at 
greater distances from the magnetic plane (see discussion 
and Figs 6-7 in Smith et al. 2006). We found that a tri- 
angular shaped occultation profile can be reproduced if 
the disk has a limited radial extent (i.e., it is ring-like) 
and/or its height perpendicular to the magnetic plane, is 
small. Another signature of a ring-like structure is if the 
ingress/egress edges of the occultation event are sharp. This 
can happen only if the disk is separated from the star. 

Additional considerations in evaluating the disk geom- 
etry in our model are first, that the ratio of secondary to 
primary line and continuum absorptions is 1.3 — 1.4. We 
will use this range as a starting guess for the ratio of the 
stellar areas eclipsed by the disk segments. The amounts 
of the absorptions in the circus computations discussed be- 
low will also translate to the individual occultation areas. 
This leads to a secondary constraint that the angles i and 
/3 must be near 90° and not much larger than 55°, respec- 
tively. Otherwise, the measured Balmer continuum absorp- 
tion cannot be attained by an occultation that is too weak 
when the disk crosses the "lower" limb of the star, as seen 
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from our vantage point. Finally, the fact that there is just a 
hint of an M-like reversion in the secondary occultation pro- 
file limits the value of i to at least than 80° ± 5°. (However, 
if the two peaks result from disk warping this argument is 
questionable.) With the constraints we imposed above as 
well as TO05's arguments, this fixes the (3 value as near 
55° ± 5°. These values are in reasonable agreement with 
estimates based on different criteria by other authors. 

For completeness, we point out that although the double 
lobes found during the secondary maximum are not treated 
explicitly in the foregoing analysis, it is clear that this disk 
segment is double-layered such that the density in or very 
near the magnetic plane is lower than in the two layers 
just outside it. We also note that for i > 80°, a period of 
1.19 days, and its vsmi of 140-162 kms"^ (e.g, Bolton et 
al. 1986, ROO, the mean radius of crOriE comes to 3.4 — 
3.9 R0. As is by now well known (e.g., Hunger, Heber, & 
Groote 1989), these radius estimates seem low for a B2 
main sequence star, but we can offer no solution to this 
apparent dilemma. 

3.2. Disk dimensions 

The dimensions of the disk segment responsible for the pri- 
mary occultation can also be inferred by appealing to Ha 
emission studies, which provide estimates of 5 — 6.4 R* for 
the outer disk dimension (Groote & Hunger 1982, Short & 
Bolton 1994, TO05), the accelerations and durations of the 
absorption lobes in Figs. 2 and 3, and the areal coverage ar- 
guments of the previous section. Our simulations also show 
that for the primary occultaion the outer radius must be 
at least 4 R* in order for the disk to have a finite thickness 
in the radial extent. Table 2 gives the dimensions of the 
absorbing disk segments in units of stellar radii as well as 
the mean densities and temperatures of our solutions. The 
radial thickness refers to the quantity Kouter - dinner- We 
note that the derived radius is consistent with that given 
from the durations of the lobe apparitions found in Figs. 2 
and 3. The errors for these parameters are about ±20%. 
For the primary absorption the results pertain to the rise 
time phase for the first half of the event, which lasts 0.04 
cycles. The shorter span for the decline (0.03 cycles) of the 
second half can best be simulated by decreasing the semi- 
height from 0.3R* to 0.2R*. These parameters also preserve 
the triangular shape of this primary occultation profile. 

The challenge for deriving the dimensions for the disk 
segment causing the secondary absorption maximum is that 
they should be large enough to ensure enough areal cover- 
age of the star and to produce the observed absorption in 
the Balmer lines and the Balmer jump. In general, lower val- 
ues of P and Kouter and larger /i values insure this. A value 
of Rinraer > 2R* is impossiblc to reconcile with the observed 
amplitude and shape of the occultation profile with phase 
in our models with the assumed /? and i angles. However, 
if one adopts the view that the extended lifetime of the oc- 
cultation is due to seeing two different radial zones of the 
same segment along a line of sight (i.e., the disk is warped), 
one can relax the distance Router up to 4R:,, (but no higher) . 
Since this constraint is consistent with the ratio of the radii 
suggested by the two acceleration rates determined in §3.2, 
we cannot rule this out. We have added this possibility in 
Table 2 by showing the admissible ranges of the inner and 
outer radii, parameter to the assumption of disk warping, 
i.e. 2R,t < Router We also point out that in the warping 



interpretation, Hmner can move in to 1.2R*, so the disk 
extends almost to the star's surface. 

According to our picture, at the middle of the first oc- 
cultation at = the whole of the disk segment is sil- 
houetted against the star. For this portion of the disk we 
may estimate the mass from the dimensions of the segment 
(area and column density) and its density. This gives a mass 
of about 5x10^^^ ^Q- If we take half the mass loss rate 
as the rate of resupply of this hemispheric disk segment 
(1.2x10-9 M0yr-i; Krticka, Kubat, & Groote 2006), the 
mass loss will replenish this disk segment in about 2 weeks. 
Because a small fraction of the wind emanating from close 
to the magnetic poles will avoid being focused into the disk 
and, further, because we may not see the full disk segment 
in azimuth around the hemisphere, the actual lifetime of 
the particles in the disk will be somewhat longerQ Note 
that the transit of a disk across the star lasts only ^10^"^ 
of this time. Therefore, leakage of disk particles either into 
the star or through its outer edge would produce a tapered 
absorption wing that is too small to be visible in our data. 

We conclude this discussion by calling attention to the 
EW minima within the general model for the disk, which 
occurred on a total of three of the seven nights of our ob- 
serving run on this star. One explanation for these min- 
ima is that a second translucent component exists layered 
over the more confined opaque disk. In this explanation, 
we would be looking through a light "haze" (with about 
1% of the column density) through most viewing aspects. 
However, little if any tail is observed in the curves of Fig. 1. 
(This fact underscores our finding above that the density 
contrast of the disk proper relative to circumstellar matter 
at moderate or high magnetic latitudes is rather high.) A 
second possible way to explain the EW minima is that they 
are "reflexes" of small incipient emissions that may still be 
detectable in the high level lines. Indeed, FEROS observa- 
tions of the Ha line at these phases (Groote 2007) exhibit 
strong emissions and an elevation of the line core at these 
phases. In this event there would be no need of a pervasive 
low-absorption component. 

3.3. Thermodynamic parameters of the disk segments 

The Balmer absorption spectra of Fig. 4 provide sufficient 
information to estimate the temperature, electron density 
and column density of the intervening disk segments. To 
make these estimates, we must assess the importance of er- 
rors in various factors used in these calculations. First, and 
importantly, the assumed areal coverages depend strongly 
on the geometrical angles derived above. Despite our agree- 
ment with earlier authors, errors in our estimates of these 
parameters propagate to the areas in complex ways and 
should be considered the largest source of our uncertainty. 
Second, it turns out that the disk microturbulence, ^, is not 
an important factor. We find almost identical solutions for 
fitting excess Balmer line and continuum absorptions for 
values of ^ in the range — 20 kms^^. We have arbitrarily 
adopted a value ^ = 10 kms~^. Third, we have used non- 
LTE model atmospheres computed from the BSTAR2006 
grid (Lanz & Hubeny 2007) and found that the photo- 



^ The fact that our column densities are averaged over the disk 
height will cause the total density and hence mass in our model 
to be underestimated. We may reasonably expect the longevity 
of particles in the disk to be a few months. 



Table 2. 



Derived absorbing disk parameters (assuming i — 80°, (3 = 55°) 
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Outer radius 


5 - 6±0.5R. 


Outer radius (2 solutions) 


2 - 4R, 


Radial thickness 


2.5 - 3i?. 


Radial thickness 


0.25 - 2.8R, 


Height 


0.3±0.07R* 


Height 


0.45±0.10R, 


Temperature 


10,500K ilOOOK 


Temperature 


12,000K ilOOOK 


Ne 


l±0.5xl0^^ cm-3 




l±0.5xl0^^ cm"^ 


Column density 


1±0.5 xlO^" cm-2 


Column density 


l±0.5xl0^'' cm-2 


Areal coverage 


0.45±0.10 


Areal coverage 


0.31±0.10 



spheric Balmer jumps for LTE and non-LTE model atmo- 
spheres differ negligibly given the adopted Tgff for cr OriE. 
This can be inferred also by examining Figure 6 of Lanz 
& Hubeny. These authors find that changes in abundances 
have little effect on photospheric temperatures in the range 
Tc = 0.01 — .003, where much of the Balmer line cores is 
formed. 

SGOl pointed out that the presence of certain low res- 
olution line aggregates in the far-UV spectra suggest the 
presence of a cool component to the disk of a Ori E having 
a temperature of about 12 750 K ± 750 K|3 In this study we 
discovered a weak but important diagnostic in the excess 
absorption spectrum of the secondary occultation (Fig. 4b), 
namely weak absorption blends at A3760 and AA3855— 60 
arising primarily from transitions of once ionized light and 
iron-group metals. Our circus syntheses show that these 
lines appear for Tdisk <12 000K over disk column densities 
of interest. 

Our primary modeling diagnostic for column density 
and areal coverage is the excess Balmer jump noted in 
Fig. 4. While potentially useful, the Balmer lines themselves 
were found to be enhanced primarily by non-LTE effects. 
Thus, whenever the excess absorption below « A3647 could 
be fit, the disk Balmer line spectrum is far too strong. We 
will therefore confine ourselves to fitting the Balmer jump 
in the following discussion. 

The electron density Ng of the disk medium can be de- 
termined by the presence of the highest-level Balmer lines 
in the occultation spectra (Fig. 4). One might first draw 
an analogy to the prominence of these lines in the spectra 
of B supergiants. For low density atmospheres, the lines 
are visible because the weak Stark wings cause less blend- 
ing. This motivated us to inquire whether the Ingiis- Teller 
(1939) criterion might apply in this instance. This condition 
gives as an estimator that the last visible line of an hydro- 
genic spectrum scales with the electron density as N~ By 
comparing the highest visible Balmer lines in the spectra 
of early-BIa stars (e.g., Bagnulo et al. 2006) and appealing 
to model atmospheres for stars of these temperatures and 
gravities (e.g., Crowther & Lennon 2006), we found that 
the visibility of H24 and H23 in occultation spectra sug- 
gest Ne values of 1x10^^ and 2x10^^ cm~^, respectively. 
However, this argument may not be precise as our spectra 
show that the strengthening of these lines occurs mainly 



^ For studies of different lines and in a different wavelength 
regime, we estimate that the uncertainty of this mean disk tem- 
perature should be at least doubled. This uncertainty does not 
propagate to the more precise temperature difference quoted 
below for opposing disk segments. 



in the line cores. Our modeling of these absorptions sug- 
gests that what is important in determining the visibility 
of the last Balmer line is the column density of high-level 
hydrogenic atoms producing the absorption. In particular, 
we find that through the Saha equation Ng is largely de- 
termined by fixing the column density at which the optical 
thickness of a high-level line is «1. For temperatures of in- 
terest (10 000-14 000 K) the density corresponding to the 
semi- visibility of H23-H24 is also = 1 x 10^^ cm"^. We 
found that the fitting of the Balmer continuum absorption 
is insensitive to the value of Ng. Therefore, we adopt this 
density for both disk segments in our circus analysis of the 
disk segment parameters. 

We adopted the densities from the above arguments to 
determine disk temperatures, column densities and frac- 
tional areal coverage factors needed from the Balmer jump 
found below A3700 in Fig. 4. To add to these constraints, 
we also invoked the condition that the areal coverage ratio 
should be 1.3 — 1.4. However, we modified this by deter- 
mining that a correction of about 10% for limb darkening 
(assuming a coefficient /i = 0.6) is necessary to account for 
the transit of the disk segments across noncentral regions 
of the star. This consideration brings us to a ratio of 1.5 for 
the fractional areas. This condition is met by the entries in 
Table 2. Our models established quickly that the hydrogen 
column density should be near 1 x 10^** cm~^. (However, the 
column lengths are not determined well enough to compare 
the radial distances of the two disk segments.) This condi- 
tion derives from the broad opacity plateau at « IxlO^"' 
cm"^ in this wavelength region for the relevant temper- 
atures. For substantially larger column densities the disk 
plasma becomes optically thick in the continuum, and the 
Balmer jump it forms begins to decrease. For smaller col- 
umn densities one runs quickly into unrealistically large 
fractional areal coverages and the Balmer jump decreases 
owing to optical thinness of the Balmer continuum. Given 
these considerations, we believe that the column densities 
we found, which are coincidentally the same for the disk seg- 
ments causing the two occultations, are accurate to within 
a factor of two. The fractional coverages we infer depend 
on this condition combined with the absorption ratio of 
1.3 — 1.4 derived above from the two occultations. 

Our estimated column density agrees well with the de- 
rived values of the electron density, 10^^ cm~'^, the inferred 
radial extent of the disk (2.5— 3 R*) for the segment caus- 
ing the primary maximum and within a factor of two of 
that for the segment producing the secondary. Note that 
SGOl found a column density of 3x10^'^ cm~^ from the 
strengths of absorption of Fe aggregates and low excita- 
tion light-metal lines in the UV. When account is taken 
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that SGOl assumed for simplicity a full occultation of the 
stellar disk, their column density is virtually the same as 
the value we derived from the assumed density and the val- 
ues given in Table 2, 3— TxlO^'' cm^^. Our solution, assum- 
ing Tcontinuum ~ 2 and disk homogeneity, likewise agrees 
perhaps fortuitously well with the peak value Tcontinuum 
»1.6 determined by TOG05 for the optical thickness viewed 
along the central magnetic plane. The latter determination 
was made for a decentered dipole that probably most ac- 
curately reflects the magnetic geometry. 

Our homogeneous disk solution overlooks the evidence 
that the segment producing the secondary occultation oc- 
curs in two semi-resolved stages. As already mentioned, one 
interpretation of this is that the a warping occurs in the in- 
ner disk, as suggested by previous studies. In any case, this 
semi-resolved structure of the occultation is noticeable in 
light curves at several wavelengths as well. It may be that 
what observers identify as "warping" is in fact due to a de- 
parture from strict coplanarity predicted for extensive Bp 
disks for which the magnetic /3 obliquity has an intermedi- 
ate value (TO05). 

Returning to the Balmer lines, we were able to evaluate 
departures from LTE in the line source functions by com- 
paring the colunui densities needed to fit the Balmer jump 
with those in the lines. As already stipulated, we assumed 
that the Balmer continuum flux is formed in the disk in 
near LTE. By comparing the column densities needed to 
fit the lines to the Balmer continua, we found for both oc- 
cultations the same approximate ratios of source to Planck 
functions: S/B = 0.006 for H9 and S/B = 0.1 for H20. 

Finally, the temperatures for the two disk segments are 
consistent with the radiative temperatures one would ex- 
pect for radiatively excited matter close to the star. There 
is no evidence in these disk segments in particular for the 
heating responsible for the presence of the variable N V res- 
onance lines (SGOl). 

4. Summary and Conclusions 

We have used high-level Balmer absorption data finely sam- 
pled around the rotation cycle to show that the distribu- 
tion of plasma in the magnetosphere of a Ori E is more 
complicated than previously believed. Most especially, the 
secondary occultation is comprised of two semi-resolved 
events. In addition, there seem to be at least two brief ap- 
pearances of absorbing matter well beyond the disk plane. 
The disk plane is dominated by two segments which oc- 
cult portions of the star some 0.4 cycles apart. We have 
no direct information from absorption profiles alone on the 
distribution of matter at other azimuths in the plane and 
so look forward to future analyses based on the Ha and/or 
H/3 emissions over the cycle. 

The occulting disk segments are sections of a ring which 
occult the star as the frozen disk corotates around the star. 
The segment causing the primary occultation lies further 
from the star (between limits of «2.5R* and 6R*) than the 
opposing segment. There are at least three arguments that 
point to the disk segments lying at different orbital radii. 
First, gas in the segment causing the primary occultation 
has a lower temperatures (10 500 K, compared to 12 000 K). 
This can be inferred first from the lack of an enhanced He I 
feature and also from the absence of weak absorptions of 
metallic-line blends referred to in §4.3. Second, the accel- 
eration of spectral migrating components across the line 



profile is higher for this segment. Third, our occultation 
analysis program indicates that one cannot reconcile the 
radii of the two segments even if one attributes the broad 

secondary absorption maximum as being due to viewing 
two occultations of a single warped disk segment. 

Perpendicular to the magnetic disk, the semiheights 

of the two segments are some 0.3— 0.45R*. It is impor- 
tant to stress that the disk of a Ori E has a smaller semi- 
height than disks for other disk-harboring Bp stars that we 
have analyzed with similar modeling techniques (e.g., SGOl, 
Smith et al. 2006). It is tempting to speculate that a disk 
should be more compressed toward the plane due to the 
pressure provided by a large magnetic field (e.g., Babel & 
Montmerle 1997). However, this cannot be the full story be- 
cause the most relevant quantity, the ratio "ry" of magnetic 
to wind energy densities is actually larger for 36Lyn than 
for a Ori E and yet the disk of the former has a somewhat 
larger height (Smith et al. 2006). Might other considera- 
tions, such as the larger radial extent of the disk (e.g. in 
the case of 36Lyn), have a bearing on this question? 

In addition to these dimensional properties, we find that 
a phase discrepancy of -1-0.018 cycles exists between the 
phase at which the central absorption transits the center 
of the hydrogen line and the phase of EW maximum. This 
agrees with a similar phase difference found by GH77 and 
calls into definition of the star's period when determined 
from different observation modes. The secondary occulta- 
tion is distinguished by the passage of two weak lobes, the 
phases of which are consistent with the semi-resolved struc- 
ture of the light curves in the literature. The interpretation 
of this structure is not clear, but it may be related to a disk 
"warping" or a corrugation of the disk plane predicted in 
the TO05 models. 

When we view the disk edge-on through the magnetic 
plane, the Balmer line absorptions are opaque, and thus the 
absorptions of the high-level members of this series increase. 
For H9 the optical depth is about 100, while it is w 10 for 
H20, and « 2 at the Balmer edge. As the disk ingresses 
and egresses there is an absence of an absorption tail in 
the equivalent width curves. This confirms the theoretical 
expectation (e.g., TO05) that the density scale height is 
very small compared to the stellar radius. This implies in 
turn that our homogeneous slab models of the disk and the 
total mass estimate of about lxl0~^° Mq are only rough 
approximations . 

How would hypothetical data of indefinitely high quality 
(signal-to-noise, spectral resolution, and cadence) improve 
on our analysis of the disk of a Ori E? We can think of the 
following applications of enhanced quality data: 

One improvement would be in the mapping of verti- 
cal density distribution and indeed the separation of the 
two density maxima associated with the secondary occul- 
tation. This would result from tracing details of the ab- 
sorption subfeatures as they migrate across the line pro- 
files. An improvement in our understanding of the vertical 
stratification should indirectly improve out understanding 
of whether the disk is warped and the inner edge of the 
secondary occultation-causing segment extends almost to 
the star or not. It will also clarify the cause of the phase 
delays that both we and GH76 have found among various 
hydrogen lines and with respect to the continuum. A reso- 
lution of this issue would allow an improved determination 
of ephemerides from fiux curves derived from different di- 
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agnostics, thereby leading to a definitive determination of 
the period and perhaps its first derivative. 

Second, the changes in the highest-level lines detected 
during the ingress and egress of the occnltations can be used 
in principle to further map the vortical density distribution 
of the disk segments because of the different optical depth 
regimes in which the various lines are formed. 

Third, highly accurate difference profiles should allow 
one to search for asymmetric absorption features in the EW 
curves that could indicate the rate at which disk particles 
leak from the disk either through the inner or outer edge. 

Fourth, increased data quality improvements (assuming 
a stable spectrograph and atmospheric transmission) can 
lead to the observation of subtle changes in disk structure. 
We suggest that these might occur on timescales of one to 
a few months or longer. Finding such changes could allow 
one ultimately to correlate them with the expected "break 
out" of disk matter through the outer edge (ud-Doula, 
Townsend. & Owocki 2006), possibly associated with X- 
ray flares from this star, (Groote & Schmitt 2004, Sanz- 
Forcada, Franciosini, & Pallavicini 2004). 
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